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Abstract 


The generalized mathematical formulation of sloshing dynamics for partially 
filled liquid of cryogenic superfluid helium II in dewar containers driven by 
both the gravity gradient and jitter accelerations applicable to scientific 
spacecraft which is eligible to carry out spinning motion and/or slew motion for 
the purpose to perform scientific observation during the normal spacecraft 
operation are investigated. An example is given with Gravity Probe-B (GP-B) 
spacecraft which is responsible for the sloshing dynamics. The jitter 
accelerations include slew motion, spinning motion, atmospheric drag on the 
spacecraft, spacecraft attitude motions arising from machinery vibrations, 
thruster firing, pointing control of spacecraft, crew motion, etc. Explicit 
mathematical expressions to cover these forces acting on the spacecraft fluid 
systems are derived. The numerical computation of sloshing dynamics have been 
based on the non- inertia frame spacecraft bound coordinate, and solve time- 
dependent, three-dimensional formulations of partial differential equations 
subject to initial and boundary conditions. The explicit mathematical 
expressions of boundary conditions to cover capillary force effect on the liquid- 
vapor interface in microgravity environments are also derived. The formulations 
of fluid moment and angular moment fluctuations in fluid profiles induced by the 
sloshing dynamics, together with fluid stress and moment fluctuations exerted on 
the spacecraft dewar containers have been derived. Results of this study have 
been widely published in the open journals. 
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I. Introduction 


For the purpose to carry out scientific experiments, some experimental 
spacecraft use cryogenic cooling for observation instrumentation and telescope, 
superconducting sensors for gyro read-out and maintain very low temperature near 
absolute zero for mechanical stability. The approaches to both cooling and 
control involve the use of superfluid liquid helium II. In this study, sloshing 
dynamics associated with spinning motions are investigated. To cover the 
spacecraft spinning motions, the Gravity Probe-B (GP-B) spacecraft has been 
chosen as the example in this study. The GP-B spacecraft adopts the cooling and 
boil -off from the cryogenic liquid helium dewar as a cryogen and propellant to 
maintain the cooling of Instrumentations, attitude control and drag- free 
operation of the spacecraft. The potential problems for cryogenic liquid in 
dewar container could be due to asymmetry in the static liquid helium 
distribution and to perturbations in the liquid-vapor interface caused by slosh 
wave excitation driven by pointing control, machinery vibration, etc. 

For the case of the GP-B spacecraft, cryogenic liquid helium II, at a 
temperature of 1.8 K, is used as the propellant. With its superconducting 
behavior, there is no temperature gradients in the liquid helium. In the absence 
of temperature gradient along the surface which drive Marangoni convection, the 
equilibrixim shape of the free surface is governed by a balance of capillary, 
centrifugal and gravitational forces. Determination of liquid-vapor interface 
profiles based on computational experiments can uncover details of the flow which 
can not be easily visualized or measured experimentally in a microgravity 
environment . 

The instability of the liquid- vapor interface surface can be induced by the 
presence of longitudinal and lateral accelerations. Slosh waves are, thus, 
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excited which produces high and low frequency oscillations in the liquid 
propellant. The sources of the residual accelerations range from the effects of 
the Earth's gravity gradient and jitter accelerations which include, atmospheric 
drag on the spacecraft, vibration of compressor, spacecraft attitude motions 
arising from machinery vibrations, thruster firings, spacecraft slew motion, 
pointing control of spacecraft, crew motion, etc. Recent study (Kamotani et 
al. , 1981) suggest that the high frequency accelerations may be unimportant in 
comparison to the residual motions caused by low frequency accelerations. 

Time -dependent dynamical behavior of partially -filled rotating fluids in 
reduced gravity environments was simulated by numerically solving the Navier 
Stokes equations subject to the Initial and the boundary conditions (Hung and 
Shyu, 1991 a,b,c; 1992 a,b,c; Hung et al., 1991 a,b,c; 1992 a,b,c). At the 
interface between the liquid and the gaseous fluids, both the kinematic surface 
boundary condition, and the interface stress conditions for components tangential 
and normal to the interface, were applied (Hung and Shyu, 1991 a,b,c; 1992 a,b,c; 
Hung et al., 1991 a,b,c; 1992 a,b,c). The initial conditions were adopted from 
the steady-state formulations developed by Hung et al (1989 a,b,c). Some of the 
steady- state formulations of interface shapes were compared with the available 
experiments carried out by Leslie (1985) in a free-falling aircraft (KC-135). 
The experiments carried out by Mason et al (1978) showed that the classical fluid 
mechanics theory is applicable for cryogenic liquid heliiim in large containers. 

In the spacecraft orbit around the Earth, the direction of azimuth angle 
of Earth toward the location of the spacecraft geometric center varies from 0° 
along the rolling axis of spacecraft to various directions in which three 
dimensional calculation shall be assijmed. 

As the spacecraft moves along the orbit, any fluid capable of motion 
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relative to the spacecraft is subject to the acceleration that arises from the 
gravity gradients of the Earth (Avduyevsky, 1984; Forward, 1982; Misner et al., 
1973) . The interaction between the particle mass of fluid and the spacecraft 
mass due to gravity gradient accelerations (Forward, 1982) are capable for the 
excitation of slosh waves and disturb the fluid system which induces the 
fluctuations of viscous stress and its moment exerted on the containers of the 
spacecraft fluid systems. In the meanwhile, the sources of residual acceleration 
of gravity jitter range from atmospheric drag on the spacecraft, background 
gravity, spacecraft attitude motions arising from machinery vibrations, 
spacecraft slew motion, thruster firings, crew motion, etc., are also capable for 
the excitation of slosh waves on the fluid containers. 

It is critically important to understand the physical and dynamical 
behavior of cryogenic helium in a rotating cylinder to effectively promote space- 
oriented missions. 

At temperatures close to absolute zero, quanttim effects begin to be of 
importance in the properties of fluids. At a temperature of 2.17®K, liquid 
helium has a A-point (a second-order phase transition); at temperatures below 
this point liquid helium (helium II) has a number of remarkable properties, the 
most important of which is superfluidity. This is the property of being able to 
flow without viscosity in narrow capillaries or gaps. 

The basis of the dynamics of helium II is the following fundamental result 
of microscopic theory. At temperatures other than zero, helium II behaves as if 
it were a mixture of two different liquids. One of these is a superfluid and 
moves with zero viscosity along a solid surface. The other is a normal viscous 
fluid. It is of great importance that no friction occurs between these two parts 
of the liquid in their relative motion, i.e., no momentum is transferred from one 
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to the other. 


It should be emphasized that regarding the liquid as a mixture of normal 
and superfluid parts is no more than a convenient description of the phenomena 
which occur in a fluid where quantum effects are important. One of these motions 
is normal and has the same properties as the motion of an ordinary viscous fluid, 
but the other is the motion of a superfluid. The two motions occur without any 
transfer of momentum from one to another. We can, in a certain sense, speak of 
the superfluid and normal parts of the fluid, but this does not mean that the 
fluid can actually be separated into two such parts. 

If fluid flow can separate helium II into the regions of the superfluid and 
the normal fluid, two temperatures zones are immediately created. A very low 
temperature zone is located at the zone of very high density concentration of the 
superfluid, while a high temperature (below 2.17® K) zone is located at the zone 
of very high density concentration of the normal fluid at the other end. The 
existence of a sharp temperature gradient at the interface between the superfluid 
and the normal fluid results in the creation of a great difference in chemical 
potential, which, in turn, induces a great reverse pressure gradient, creating 
the environment of isothermal fluid distribution everywhere throughout the 
cylinder and homogenous distribution of superfluid density concentration. This 
illustration of the possible separation of superfluid from normal fluid of helium 
II means that there is in reality, no way for anyone to achieve the separation 
of the superfluid from the normal fluid of helium II. In other words, in 
considering the dynamical behavior of helium II In a large rotating cylinder, a 
mixture of the superfluid and the normal fluid without separation of the two 
fluids is accounted for in the model computation. Density concentration of 
superfluid is a function of temperature, which is also true for the surface 
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tension and viscous coefficient for helium II (Wilks, 1967; Hoare et al., 1961; 
Hung, 1990; Hung and Lee, 1992). In fact, the experiments carried out by Mason 
et al. (1978) showed that the classical fluid mechanics theory is applicable for 
cryogenic liquid helium in a large container. In this study, the theory of 
viscous Newtonian fluids is employed with modification of transport coefficients 
adjusted by normal and superfluid density concentration which is a function of 
temperature . 

II. Functions of Scientific Observation and Spacecraft Motions 

The GP-B spacecraft is a sun synchronous Earth satellite orbiting at 650 
km altitude directly over the poles. The functions of scientific observation for 
the GP-B spacecraft and its motions are illustrated as follows: 

The GP-B spacecraft is a relativity gyroscope experiment to test two 
extraordinary, universified predictions of Albert Einstein's general theory of 
relativity (Wilkinson et al., 1986; Stanford Relativity Gyroscope Experiment, 
1986), By using gyroscopes (ones with electrically supported spheres, spinning 
in a vacuum and others which utilize the spins of atomic nuclei, circulating 
sound waves, and even circulating laser beams), the GP-B measures two distinct 
space- time processes, frame dragging and the geodetic effect , which gradually 
changes its directions of spin. In these gyroscopes, the underlying principle 
is that rotating systems, free from disturbing forces, should remain pointing in 
the same direction in space. 

The physical meaning of "the same direction in space" can be explained from 
the differences in the Newtonian and Einsteinian concepts of the universe. For 
Newton, the answer was easy, in that space and time were absolutes. A perfect 
gyroscope set spinning and pointed at a star would stay aligned forever. Not so 
for Einstein, for whom space and time is warped. A gyroscope orbiting the Earth 
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finds two distinct space -time processes, frame drajzping (rotation of space -time) 
and the geodetic effect (curvature of space-time), which will be measured by the 
GP-B. 

As to frame -dragging for disturbances in the rotation of space- time, 
calculations based on a gyroscope in polar orbit at 650 km altitude should turn 
with the Earth through an angle amounting, after one year, to 42 milliarc- seconds 
(1 milliarc -second - 10"^ arc -second; 3.6 x 10^ arc -seconds - 1®, angular degree) 
(Everitt, 1990). As to the geodetic effect for disturbances in the curvature of 
space -time, the gyroscope's motion for rotation in the orbit plane at a 650 km 
altitude, from relativistic correction to the complicated motion of the Earth- 
Moon system around the Sun, should turn 6,600 milliarc-seconds per year (Everitt, 
1990) . 

Furthermore, one of the GP-B's most important tasks is to investigate the 
gravitational action of moving matter in the universe. In other words, matter 
moving through space- time can be thought of as creating a new force- 
gravitomagnetism, which is similar to electromagnetic force created by an 
electric charge moving through magnetic fields (Misner et al., 1973; Everitt, 
1990). The frame -dragging measurement detects this force and fixes its scale. 

In order to insure these extremely precise and accurate measurements and 
gyroscope operations, near zero temperature is required for mechanical stability 
of the instrtunent, preservation of the lead bag magnetic shield, shielding the 
gyroscopes against nongravitation disturbances , and for reading their direction 
of spin. Near zero temperature cryogenic liquid helium II (1.8®K) has been 
chosen to serve this purpose (Mason et al,, 1978; Hung, 1990), 

The GP-B instrument comprises four gyroscopes and a reference telescope 
sighted on Rigel, a bright star in Orion, In polar orbit, with gyro spin 
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directions also pointing toward Rigel, the frame dragging and geodetic effects 
come out at right angles, each gyroscope measuring both. The GP-B attitude 
controller with a fine pointing system inside the dewar probe will constantly 
assure the operation of the spacecraft within the range of the scope of each 
experiment through the proper control of thrust. 

Thrust comes freely from the dewar 's boil-off helium vapor, in contrast to 
a traditional gas jet control system which uses on-off valves to save gas. Here, 
the copious supply of vapor allows the use of proportional thrusters : pairs of 
opposed nozzles adjusted continuously for much smoother control. Six thruster- 
pairs, suitably arrange, cover all three control functions with redundancy in 
case of failure. To take advantage of the characteristics of viscous-free 
superfluid helium II under X transition point, a porous plug is used to control 
the helium boil-off and to produce the fountain pressure to serve as their thrust 
jet for the spacecraft propulsion system. 

To comprise these functions of scientific observation, GP-B stores its 
gyroscopes, telescope, probe mass and others in the center core of the dewar 
probe surrounded by the cryogenic helium II liquid. The dewar container of the 
GP-B is spinning with a rotating rate of 0.1 rpm during normal operation. As the 
telescope is constantly sighted on Rigel, gyro spin directions also pointing 
toward Rigel during the spacecraft moving around the polar orbit. The GP-B with 
its rotating dewar, there is no slew motion involved in the spacecraft motion. 

III. Basic Charcteristics of Gravity Gradient 
and Gravity Jitter Accelerations 

Any fluid element inside the on-orbit spacecraft fluid system is subject 
to the acceleration that arises from the gravity gradient of the Earth 
(Avduyevsky, 1984; Forward, 1982; Misner et al., 1973; Hung and Pan, 1993; Hung 
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et al. , 1993 a,b,c). Once the spacecraft orbit is fixed, the orbit period is 
determined and the basic structure of the gravity gradient acceleration also can 
be calculated. However, gravity gradient acceleration acting on each fluid 
element inside the on-orbit spacecraft fluid system is different dependent upon 
the distance of the location of the fluid element to the geometrical center of 
the spacecraft and its direction toward the location of the center of the Earth. 
This acceleration can only be calculated based on the non- inertia frame of 
spacecraft bound coordinate. Thus, the coordinate system shall be transformed 
from ordinary inertia frame coordinate to non- inertia coordinate. 

(A) Orbit Motion of Spacecraft 

Let us consider the case of the GP-B spacecraft, which is the Earth 
satellite orbiting at 650 km altitude directly over the poles, the orbit period, 
Tq can be computed from following expression: 


x^=2« 


R. 


3/2 


RgST^^ 


(3-1) 


where denotes the radius of Earth (- 6373 km); R(., the radius of the circular 
orbit (- Rg + h - 7023 km); h, orbit altitude (- 650 km); and g^, Earth gravity 
acceleration (- 9.81 ™/s^) - For the case of the GP-B spacecraft, the orbit 
period - 97.6 min, and orbit rate n - ^/Tq " 1.07 x 10“^ 

As the spacecraft is orbiting around the Earth, the azimuth angle of the 
Earth, toward the location of the spacecraft geometric center varies with 
respect to time. At time t - 0, the rolling axis of the spacecraft is aligned 
with the radial direction of the Earth's center to the spacecraft geometric 
center. Assuming the spacecraft rolling axis is linearly turning arotind 0® to 
360® in the orbit period, of the spacecraft when the spacecraft is orbiting 
around the Earth, This is particularly true for the case of the GP-B spacecraft. 
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Without the spacecraft slew motion, the azimuth angle (^o) defined as 



where is the spacecraft orbit period [defined in Equation (3-1)]; and t is the 
time measured from the instant when the direction of the spacecraft rolling axis 
is aligned with the radial direction of the spacecraft geometric center to the 
center of the Earth. 

(B) Spinning and Slew Motions of Spacecraft 

For the purpose to carry out wide-range observations, some scientific 
spacecraft requires slew motion with respect to the mass center of the 
spacecraft. For the case of the spacecraft slew motion, azimuth angle, shown in 
Equation (3-2), shall be modified through the coordinate transformation of slew 
motion when the spacecraft is orbiting around the Earth. 

Let us assume that the slew motion starts with the center located at the 
mass center of the spacecraft. Let us choose cartesian coordinate (x" , y" , z") 
with z"-axis along the axis of the dewar container (see Figure 1). At time t - 
0, the radial vector from the center of the spacecraft to the center of the 
Earth lies on the x^-z" plane of the cartesian coordinate chosen. The azimuth 
angle ^ is defined as the angle between the radial vector tp and the z"-axis. 
Rotation matrices for spinning and/or slew motions along the x"-, y"- and z"-axes 
can be expressed as 


10 o' 


cosQyt 0 -sinco^t 


cosu^t sinw^t 0 

0 cosco^t sinOj^t 


0 10 


-sino)^t cosQ^t 0 

0 -sinUj^t cosOjjt 


sincjyt 0 cos(Oyt 

» 

0 0 1. 


respectively. Here, cjy and denote angular velocity of slew and/or 

spinning motions along the x"-, y"- and z"-axes, respectively. Radial vector ip 
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in cartesian coordinate without slew and spinning motion is 


0, -cosil^^J (3-3) 

With an execution of spinning motion along the z"-axis only, radial vector r^. 
becomes 




coso^t sino^t 0 
-sinw^c cosQ^t 0 
0 0 1 


sinilfa, 

0 


-COSllTg^ 


» [sini^g^coscij^t, -simlTg^sinu^c, -cosilTg^] (3-4) 

With an execution of slew motion along the y"-axis only, radial vector r^ becomes 


cosQyt 0 -sinci)j,t 


’ sinijr^^ 

0 10 


0 

sincOyt 0 cosOyt 


-cosilr^^, 


- [sin(Tlr^o + «yO , 0, -cos(i|rgj, + (ayt) ] (3-5) 

With an operation of slew motion along the x"-axis only, radial vector r^ becomes 


■c-x 


0 cosOjjt sinOjjt 
0 -sinOjjt cosojjjt 


sinijfg^ 

0 

|-cosi|rg„j 


= [sini|;^, -cosi|rg^sin<Oj^c, -cosTlr^oCosQ^t] (3-6) 

In other words , radial vector will be modified from the mathematical 
expression shown in Equation (3-3) to (3-4), (3-5) and (3-6) for the slew and/or 
spinning motions along the z"-, y"-, and x"-axes alone, respectively. In 
particular, for the case of slew motion along the y”-axis, comparison between 
Equations (3-3) and (3-5), it shows that the azimuth angle will be modified as 
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= ’►so “yt 


(3-7) 


For the successive operations of the spacecraft from spinning motion along 
the z"-axis, then slew motion along the y"-axis, and then slew motion along the 
x"-axis, radial vector fe results 


■c-x.y.x 


1 

0 

0 


0 1 
cosQj^t sinOjft 
-sinQjjt coscjjjt 


cosUyt 0 
0 1 
sinWyt 0 


-sinWyt 

0 

COS(i>j,t 


cosco^t 

• -sincOjjt 

0 


SixiQ^t 0 
COS<i>^ 0 
0 1 


’ sinilTg^, 
0 


[-cosi|r^j 


(3-8) 


In addition to the modification of the azimuth angle made by the spacecraft 
slew motion through the formulation of coordinate transformation, shown in 
Equation (3-3) to (3-8), accelerations are also induced to activate on the fluid 
mass in the dewar container. Accelerations acting on the fluid particle in the 
dewar induced by the slew motion of the spacecraft with the coordinate fixed at 
the spacecraft center of the mass is as follows (see Figure 1) : 

•« 

ifp = (^xiQxRp) + d,xRp + 2S>x^ (3-9) 

where ^ denotes the position vector of the fluid particle in the dewar container 
relative to the body frame of the spacecraft; angular velocity of the 
spacecraft body frame; a, angular acceleration of the spacecraft body frame; and 
velocity of the fluid particle relative to the spacecraft body frame. 

As we indicated earlier, let us assume that the slew motion starts with the 
center located at the spacecraft mass center, cartesian coordinate (x' ' , y' ' » 
z'') is chosen with origin located at the spacecraft mass center. Let us also 
assvime that x''-z'' plane intersects the center of Earth and the spacecraft mass 
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center. In other words, azimuth angle of Earth toward the spacecraft mass center 

lies in the x''-z'' plane. Slew motion is along both the x' ' - andy''- 

«« 

coordinates. Thus, 0) and ®syf 0) » ^ <^ne to slew motion 

becomes 


•f 

Ry 


sy^z 

^p,slew * 

Ryll 

R^u 

S 


slew 



(C) Coupling for the Accelerations of Spinning and Slew Motion of Spacecraft 
For some particular reasons required in the scientific spacecraft, it might 
be faced with the situation that both spinning and slew motions are needed 
simultaneously. To encounter this case, the following formulations are made to 
deal with coupling for the accelerations of spinning and slew motion of the 
spacecraft: 





slew end spinning 


gyRjf) — Ctj3pR^~2Cl)^V^ 

+2(o^Vj,-o> 


sy^x> 


J slew 


+ 


{“sy^*-«z^y) “r + “r^x + 2(0*V;* 


'spinning and coupling 


( 3 - 11 ) 


where and denote angular velocity and angular acceleration, respectively, 
of spacecraft spinning motion along the z-axis. 

For the case of the GP-B spacecraft, there is no slew motion and the 
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spinning is the only acceleration acting on the spacecraft fluid system. 
Acceleration due to spacecraft spinning motion becomes 


spinning 





n 

3 



spinning 

0 


(3-12) 


i spinning 


To convert the expression of Equation (3-12) in cartesian coordinate to 
cylindrical coordinate, by using the relationships of (R^, Ry) - (rcosd, rsind) 
and (Vjj, Vy) - (Uj. cosO - sintf, u^ sintf -h u^ costf), Equation (3-12) becomes 
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-rcos0ti),^-r3in0c!)^-2 (u^sin0 + ujcos0) u. 

R — 

^p, spinning 


spinning 

-rsin0(i)j^+rcos06^ + 2 (u,cos0-U0sin0) 

0 


(3-13) 

spinning 


and 


^p. Spinning * 

Rr 

3 

i?j^cos0 +i?ySin0 
-i?j^sin0 + i?yCO30 

- 

-ra>,2-2u9&)^ 
r(b^ + 2Ur&>^ 

(3-14) 



Spinning 


spinning 

0 

Spinning 


Accelerations induced by spacecraft spinning motion alone becomes 


■^r 


Rr 


ro>/ + 2ii9W^ 


3 — 


3 

-r6^-2u^oj. 


\ spinning 

Rz 

spinning 

0 


(3-15) 


Spinning 


(D) Gravity Gradient Acceleration 

The gravity gradient acceleration acting on the fluid mass of spacecraft 
can be shown as 


[3 f^-d] (3-16) 

where a^g denotes gravity gradient acceleration vector; the vector (not a unit 
vector) from the fluid element to the spacecraft geometric center; ^ unit 
vector from the spacecraft geometric center to the center of the Earth; and n, 
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the orbit rate (see Figure 1). 

It is assumed that the gravity gradient exerted on the geometrical center 
of the spacecraft orbiting around the Earth on its specified orbit is zero. In 
other words, all the gravity acceleration exerted on the spacecraft is nothing 
but the gravity gradient acceleration which is defined in Equation (3-16). In 
this study, we are interested in investigating how gravity gradient acceleration 
affects the dynamical behaviors of cryogenic fluid elements of helium. 

For the convenience of mathematical calculation, let us describe all the 
parameters involved in Equation (3-16) in terms of cartesian coordinates. In 
order to match with the computer simulation, mathematical derivation are 
considered in the first quadrant. Figure 1 illustrates the geometrical 
relationship of the parameters shown in Equation (3-16). 

Let us consider the fluid element of interests, m, located at (r, 5, z) in 
cylindrical coordinates and at (x, y, z) in cartesian coordinates. The origin 
of the two coordinate systems is located at the center bottom of the dewar tank. 
The slew and/or spinning motions, mentioned earlier, are executed at the 
spacecraft mass center with cartesian coordinate (x" , y" , z"). The geometry 
center or the spacecraft mass center is located at z - Lg- As |d| (not an unit 
vector) is much smaller than the distance between the location of the GP-B 
spacecraft geometric center to the center of the Earth, r^ (an unit vector) 
through the GP-B geometric center and r^ (an unit vector) through the fluid 
element, m, is basically the same. Assume that vector r^ lies in the x-z plane 
of the cartesian coordinate. 

Radial vector r^ with the modification of slew and/or spinning motions 
along the x"-, y"-, z"-axes have been derived in Equations (3-3) to (3-8). Based 
on the relationship between coordinates (x, y, z) and (x**, y", z") 
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(3-17) 


X 


i 0 0 


x" 

y 

= 

0 10 


y" 



O 

O 

M 


z" 


Vector in (x, y, z) coordinate becomes 


d= [-rcos0, -rsinO, -(z-Lj,) ] (3-18) 

Substituting Equations (3-4) and (3-18) in (3-16), non-inertia frame 
expression of gravity gradient acceleration with spinning motion in z-axis 
becomes 






’3 [-rslnf cos^ J sin^j^cosw^t^rcosd 

-3 [>rBini|rj^cos (6-t-tt.C) ♦ (r-L^) 8in^^sinca,t>rsin6 

-3 [-rsinfgj^co8(6^»,t> > (z-r^) cost (z-L^) 


(3-19) 


Substituting Equations (3-5) and (3-18) in (3-16), non- inertia frame 
expressions of gravity gradient acceleration with slew motion in y-axis becomes 




3 [ -X3in»|fcos0 + cosi|r (z-Lg) ] sinijr + rcos0 

^99tY 

= xi^ 

rsin0 

^gg, z. 

sl9w in y-axis 

-3 [ -rsini|rcos0 +cost|r (z-L,,) ] cosi|r + (z-L^) 


( 3 - 20 ) 


where ^ + Wyt. 

Substituting Equations (3-6) and (3-18) in (3-16), non-inertia frame 
expressions of gravity gradient acceleration with slew motion in x-axis becomes 






n» 


3(-rco808in^j^<K;o8^^[r8in6>,t8lnd^co8G>j^C (z-LJ ]}8in^^>rco8d 
-3{-rco868in^j„+co8^ga) J8ino,t8in6+C08<iijjt(z-L^,) ] /cosijrj^sinUjrZ^rzlnO 
-3(-rco8d8inf j^-k:o8^^ (z-L^) ])co8f^coa<i)j,£:+ (z-Lg) 


( 3 - 21 ) 


The relationship for the coordinate transformation from cartesian to 


cylindrical coordinates for any vector F (such as velocity or force vectors ) in 
non- inertia frame of spacecraft bound coordinate can be shown as 


(3-22) 


f^r] 


COS0 sin0 0 


\P.] 


= 

-sin0 COS0 0 


Py 



0 0 1. 




Thus, the gravity gradient acceleration located at (r, B , z) can be computed from 
that located at (x, y, z) , shown in Equations (3-18) to (3-21), from the 
following relation: 



'a 


COS0 sin0 0 


a 


^779 r 
^sry-8 

= 

-sin0 COS0 0 


97 

^79. y 




O 

o 


^77. X. 


(E) Jitter Accelerations 

In addition to gravity gradient acceleration acting on the fluid element 
of on-orbit spacecraft fluid systems, there is another acceleration of gravity 
jitter also exerted forces on the fluid systems. The sources of residual 
acceleration of gravity jitter range from slew motion of spacecraft, atmospheric 
drag on the spacecraft, background gravity, spacecraft attitude motions arising 
from machinery vibrations, thruster firings, crew motion, etc. , are also capable 
for the excitation of slosh waves in spacecraft fluid systems (Kamotani et al., 
1981; Hung and Shyu, 1991 a,b,c; 1992 a,b,c; 1993 a,b,c; Hung et al., 1992 
a,b,c) . 

Among all of the varieties of jitter accelerations listed, accelerations 
induced by slew motion of the spacecraft dominate the forces activated on the 
spacecraft fluid systems. Two coordinate systems (cylindrical and cartesian) 
chosen in this study are (r, z) with corresponding velocity components (u^, 
^99 for cylindrical, and (x, y, z) with corresponding velocity components 
(u^, Uy, u^) for cartesian coordinates. The origin of these two coordinates are 
located at the central bottom of the dewar tank, as shown in Figure 1. The 


16 



spacecraft center of mass, or the geometric center of the spacecraft is located 
at (x^,, y^,, z<j) - (0, 0, L^) . The relationships of the coordinate, velocity and 
the force between cartesian and cylindrical coordinates are 
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(3-26) 


In the derivation of acceleration induced by the slew motion of spacecraft, 
the coordinate system (x**, y”, z”) is fixed at the spacecraft center of the mass. 
The relationships of the coordinate, velocity and acceleration between 
expressions with the origin located at the spacecraft center of the mass (x" , y" , 
z") and origin located at the center bottom of the dewar tank (x, y, z) are 
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(3-30) 
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A detailed expression of [^, Ry, shown in Eqxiation (3-10) of 
this report. Jitter acceleration is a simunation of acceleration Induced by slew 
motion and others, such as atmospheric drag on the spacecraft, spacecraft 
attitude motions arising from machinery vibration, thruster firing, crew motion, 
etc. Thus, jitter acceleration can be expressed as 
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where f is the jitter frequency (Hz) imposed on the fluid systems of the 
spacecraft. 

IV. Non- Inertia Frame Mathematical Formulation of Fundamental Equations 
Dynamical behavior of fluid elements inside the on-orbit spacecraft fluid 
systems are strongly modified by the gravity gradient and gravity jitter 
accelerations . In order to accommodate the impact of gravity gradient 
acceleration, in particular, on the on-orbit fluid motion, one has to consider 
non- inertia frame of the spacecraft bound coordinate rather than adopting inertia 
frame coordinate used in ordinary fluid mechanics formulation. 

Consider a closed circular cylindrical dewar of radius, a, with height, L, 
which is partially filled with cryogenic liquid helium, and rest of the ullage 
is filled with a helium vapor. Angular velocity of rotating cylinder is w. 
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Density and viscosity of liquid helium and helium vapor are Pl, Pv, and /iy, 
respectively (Mason, 1978). Let us use cylindrical coordinates (r, z), with 
corresponding velocity components (u, v, w) , and corresponding residual gravity 
acceleration, such as gravity gradient components gravity 
jitter components derivation of the governing 
equations, accelerations induced by the spinning motion of the spacecraft is 
included in the formulation. The rest of the acceleration such as slew motion, 
atmospheric drag on the spacecraft, spacecraft attitude motions arising from 
machinery vibrations, thruster firing and others, are included in the jitter 
acceleration, shown in Equation (3-31). The governing equations for non- inertia 
frame of spacecraft bound coordinates can be shown as follows: 

(A) Continuity Equation 
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(B) Momentum Equations 
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(4-5) 


V==A4-(r^) ^ 


I dr dr ^02 dz^ 

In these formulations, 2wjV and ZcjjjU are the Coriolis acceleration, rw^^ is 
the centrifugal acceleration, and rci^ is the angular acceleration induced by the 
spinning motion of the spacecraft. 

In the computation of fluid forces, moment, viscous stress and angular 
momentum acting on the container wall of the spacecraft, one has to consider 
those forces and moment in the inertia frame rather than the non- inertia frame, 
in particular for the case of the spinning motion in the z-axis. To show an 
example, one has to transform those vectors from the non- inertia frame to the 
inertia frame for the case of spinning motion in the z-axis. 
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where prime symbol denotes vectors in the inertia frame while those parameters 
without the prime symbol indicate vectors in the non- inertia frame. 

V. Initial and Boundary Conditions of 
Spacecraft Fluid System in Microgravity Environment 
Governing equations of the fluid motion in on-orbit spacecraft fluid 
systems for non- inertia frame of spacecraft bound coordinates have been 
illustrated in Equations (4-1) to (4-6). These equations shall be combined with 
the characteristics of gravity gradient and gravity jitter accelerations as 
formulated in Equations (3-1) to (3-31). Initial and boundary conditions shall 
be introduced to accommodate solving fluid motion in on-orbit spacecraft fluid 
system for non- inertia frame coordinate (Hung et al, , 1990 a,b,c; 1991 

a,b,c,d,e,f,g,h,i,j). 
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Let the profile of the interface between gaseous and liquid fluids be given 

by: 

nit, r, z) - 0 . (5-1) 

The initial condition of the profile of the interface between gaseous and 
liquid fluids at t - t^ is assigned explicitly, and is given by: 

nit - to, r, z) - 0 (5-2) 

A set of boundary conditions has to be supplied for solving the equations. 
These initial interface profiles used in this study have been given explicitly 
through the steady state computations made by Hung and Leslie (1988) and Hung et 
al (1989 a,b,c,d) which were checked by the experiments carried out by Leslie 
(1985). These boundary conditions are as follows: 

(1) Along the container wall, the following three boundary conditions 
apply: 

(a) Interface between solid and liquid fluid: No -penetration and no-slip 

conditions assure that both the tangential and the normal components 
of the liquid velocity along the solid walls will vanish. 

(b) Interface between solid and gaseous (vapor) fluid: Similar no- 

penetration and no-slip conditions as that shown for interface between 
solid and liquid fluid will apply. 

(c) At the location of solid- liquid-gaseous (vapor) three phases 

interface: No-penetration, but not no-slip condition apply. This 

will assure that normal components of liquid and vapor velocities 
along the solid wall vanish, and allow a slipping flow of liquid and 
vapor fluids along the solid wall at three phase interface location. 
The velocity of slipping flow at this location is governed by the 
adhesive forces between fluids (liquid and gaseous) and solid walls. 
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Also, ac this location of three phase interface, a constant contact 
angle is present in which the behaviors of wet or dry contacts are 
detennined by Coulomb interaction between the fluids (liquid and 
vapor) and the surface phenomena (material and roughness) of solid 
walls . 

(2) Along the interface between the liquid and gaseous 

fluids, the following two conditions apply: 

(a) Kinematic surface boundary condition: The liquid (or gaseous) surface 
moves with the liquid (or gas) which implies 

012 n ( t=c,, r,d,z) 

(b) Interface stress condition: Across the liquid-vapor interface, the 
stress must be continuous. Based on Landu and Lifshitz (1959), the 
stress across the liquid-vapor interface can be expressed as 

(Ps-F^)12^-[(x^j) + -i)i3^ (5-4) 

where and R 2 are the radius of curvatures of two major axes at the point of 
interests on the surface of the liquid-vapor interface. 

The expressions of radius of curvatures R^ and R 2 in cylindrical 
coordinates from differential geometry can be shown as 

1 ^ 1 . EN-2FM^GL 

R 2 EG-F^ ^ ’ 

where the relationship of cartesian and cylindrical coordinates for the curved 
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Here, the configuration of the liquid-vapor interface is z - H (t - t<,, r, S); 
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[xg, yg, Zg\ = [x, y, z] 
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Simplifying Equations (5-7) to (5-11), one obtains 


EG-F^ = r^{l + H/ + — 


E = 1 H/, F = G ^ 
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D rD rD 

Substituting these relations to Equation (5-4), the radius of the curvature 
on the curved surface of the configuration of liquid-vapor interface in 
cylindrical coordinates can be expressed as follows: 


— + — = (r~) + — ( — )] 


( 5 - 15 ) 


Here, in Equations (5-4) and (5-15) 


auj 3u. 2 3u*. 


is the viscous stress tensor; ft, the viscous coefficient of the first kind; f, 
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the viscous coefficient of the second kind; P, the pressure; dr, the surface 
tension of the liquid-vapor interface; and n^, the unit vector normal to the 


interface; and the Dirac delta function. Also, subscripts G and L denote 
conditions at gaseous and liquids fluids, respectively, across the liquid-vapor 
interface . 

The fluid stresses across the liquid-vapor interface can be decomposed to 
the components normal (n^, a unit vector) and tangential (t^, a unit vector) to 
the interface. For the component tangential to the interface, one can take a dot 
product of a unit vector tangential to the interface, t^, to Equations (5-4) and 
(5-15), which leads to 

(5-16) 

since n^^t^ - 0. 

For the component normal to the interface, one can also take a dot product 
of a unit vector normal to the interface, n^, to Equations (5-4) and (5-15), 
which leads to 


Pg-Pl- ( 5 - 17 ) 

For components normal to the interface along the (r, z) directions in 
cylindrical coordinates can be obtained by taking dot products of n^, n^, 
separately to Equations (5-4) and (5-15), which are expressed as 


- Pl) 
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where (n^, n^, n^) is the unit vector normal to the interface in cylindrical 
coordinates (r, 9 , z) , 
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r ‘ dz 




(5-18) 


For a special case of axial synunecry, component normal to the interface, 
shown in Equation (5-17), can be simplified and becomes 


Pa - ^L - 


£i— 1 (5-19) 

r dr 1 

(1 + 4 )") " 

since 3/39 - 0, - 3H/3r - dz/dr - - 0 and D • (1 + for the case 

of aucial symmetry. 

VI, Characteristics of Slosh Wave Induced Fluctuations 
in Fluid Moment and Angular Momentum 

Slosh wave induced fluctuations in the fluid system of the rotating dewar 
introduce time -dependent disturbances in moment and angular momentum of 
spacecraft fluid system. In this study, there are induced angular velocities 
along the yawing, pitching and rolling axes due to the fluid motion inside 
rotating container. These angular velocities in yawing, pitching and rolling 
axes, caused by the fluid flows in a partially liquid- filled container, readjust 
the angular velocity in rolling axis. 

In order to accommodate the spacecraft dynamics of yawing, pitching and 
rolling, cylindrical coordinates (shown in Figure 1) of rotating container is 
transformed into cartesian coordinates based on (x, y, z) - (rcos3, rsiniS, z) 
with corresponding velocity components (V^^, Vy, V^) - (ucostf - vsintf, usin^ + 
vcos3, w) . If spacecraft is rotated with respect to mass center at (r^., ^c) 

in cylindrical coordinates, location of mass center in cartesian coordinates 
becomes (x^, y^, z^) - (rcCOstfg, r^sinS^, Zg) . Induced angular velocities Uy, 
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in cartesian coordinates becomes 
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As the velocity components are given by 
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the relationship between the components of induced angular velocity and flow 
velocity in cylindrical coordinates can be expressed in the following 
formulation: 
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^xx [ -(z-z^) sin9 

■^xy = - iz-z^) cosd [{rsinB-r^inQJ)^+ (,z-z^)^]~^ 
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i 0 i 

For the case of the GP-B, the axis of rotation is always fixed at the point 
of proof mass which is located at the geometric center of the dewar at (Xg, 7cfZc 
— (0, 0, Lc) where - 1/2 L and L is the height of the dewar (see Figure 1) . 
By using the computed results of induced angular velocity shown in Equation (6- 
3), one can compute the angular momentum (H^^, Hy, H^) as follows: 
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The moment of spacecraft can be computed from the time rate of change of 
the angular momentvun , i . e . , 
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(6-5) 


where cjy, a?^) denotes pitching, yawing and rolling angular velocities 

of spacecraft in inertia frame. 

VII, Mathematical Formulation of Fluid Stresses 
and Moment Fluctuations Due to Slosh Waves 
For the purpose of considering large amplitude slosh wave activated fluid 
stresses exerted on the solid walls of the dewar, the fluid stresses are 
decomposed into the tangential and normal components to the walls which can be 
expresses as follows: 


_ / du. dug \ ^ ^ 




(7-1) 




^[dx^ axj 




(7-2) 


where denotes the tangential component of fluid stresses; 11^, the normal 
component of fluid stresses; P, the thermodynamic pressure; u^, fluid velocity 
in a direction; t^,, unit vector tangential to the wall; n^, unit vector normal 
to the wall; /x, the molecular viscosity coefficient of fluid; and the Dirac 
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delta function. Subscripts a and /9 imply the directions of flow fields. 


Figures 2(A) and 2(B) show the geometry of GP-B dewar propellant tank in 
both r-z and r-^ planes, respectively. In order to make the computation of fluid 
stresses match the geometry of the dewar tank, mathematical formulations have 
been divided into three sections: (A) Top wall (dome) section, (B) Bottom wall 
(dome) section, (C) Cylindrical section (including probe section of the inner 
wall of dewar) and (D) Baffle plates section. There are several plates of baffle 
inserted in the dewar [see Figure 3(A) and 3(B) for baffle board illustration]. 
Baffles with the shape of hollow circular plate with an inner radius R^, and an 
outer radius R 2 are installed along the probe column of the dewar located at z 
- L^, where i - 1, 2, . , , , n plates and the thickness of each plate is d. Figure 
4 shows the GP-B dewar container equipped with probe and baffle boards. 

(A) Top Wall (Dome) Section: 
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(B) Bottom Wall (Dome) Section: 
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(C) Cylindrical Section (including probe section of the inner wall of dewar) : 
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(D) Baffle Plates Section: 

*-|7> 

-If) a-13) 
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where (f> is the azimuth angle of the dome; and and denote tangential 
stresses in r-z and r-5 planes of the dewar, respectively. Velocity components 
in cylindrical coordinates of (r, S, z) are shown as (u, v, w) . 

The stress distribution shown in Eqs. (7-3) to (7-14) can be integrated 
with respect to area and obtain the tangential and normal forces on top wall, 
bottom wall, and cylindrical sections of the dewar. 

(A) Top Wall (Dome) Sections: 
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(B) Bottom Wall (Dome) Section: 
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(C) Cylindrical Section (including probe section of the inner wall of dewar) 
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(D) Baffle Plate Section: 
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where and denote the radius and the azimuth angle of the dome section, 
respectively; and and Rg 2 are the radii of the inner and outer walls of the 
cylindrical section of the dewar (see Figs. 2 and 3). 

In order to accommodate the spacecraft dynamics of pitching, yawing and 
rolling, cylindrical coordinates (shown in Figs. 2 and 3) of the rotating 
container is transformed into cartesian coordinates based on (x, y, z) - (rcos^, 
rsintf, z) with corresponding velocity components (v^, Vy, v,) - (ucos5 - vsind , 
usin^ + vcosd , w) . For the case of the GP-B Spacecraft, the axis of rotation is 
always fixed at the point of proof mass which is located at the geometrical 
center of the dewar at (Xc y^, z^) - (0, 0, L/2) , where L is the height of the 
dewar (see Figs. 2 and 3). To fulfill this goal, stress distributions, shown in 
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Eqs. (7-15) to (7-18), have to be recalculated in the (x, y, z) directions (see 
Figs. 2 and 3). 
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The moment of stress force acting on the dewar wall of the container can 
be computed from the cross product of moment arm, [which is the perpendicular 
distance from the location of proof mass of the spacecraft to the total forces] , 
and the total forces. Components of the moment can be computed from the 
following formulations : 
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+U I (He) fftiij+// (He) Attorn Wiji]j2dcos2<|>-cos*0d<|>d0 

-{/ / (n^) TOP (n„) BotCo« „.ji]i?dCOS^(|)-Sin0-COS0d(t>dO 

*lff (n,)f^„cyiipd«cos0+|/ (n„)j^„cyii«d«'sin0]iecicos0d0cfz 
{ + [ JJ (H,) 5 ;;®t«cyii^d«COS 0 H- |/(n„)o„cercyiind«Sin 0 ]i?| 2 COS 0 d 0 dz 

-Uff- (n^) Top wall*// (He) B~oltom Will •J?jcos2(j>*sin<j>-cos0-sin0d(j)d9 
- (n^) it Wii -/ f (He) Ittoa, „au]Rlc03^-3±n^ed<^dQ 

■{[ / (n„) rop Wii^ <nj ^,,,]i?^cos 3 <l)-cos 0 -sin 0 d(t>d 0 

J/m«rCyiindar^^^® 

Outer cylinder^ 1110 + //®.> outer cylind0z^°^^ ] -Re2Sin0d0dz} (7-24) 

where denotes the height of the dome (see Fig. 2 ). 

After the integration of Eqs. ( 7 - 19 ) thru ( 7 - 24 ), components of the moment 
can be rewritten as the following expression: 


[MJ lLyF,-U,-L,)Fy] 
My = U,-L^) 

Mz\ L^y-LyF^ J 


( 7 - 25 ) 
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where L^, Ly» and L* denote the components of the moment arm along x, y, and z 
axes , respectively . 

By using the relation of - 0, or Ljj + Fy Ly + (L^ - L^) • 0, 
moment arms of the moment of fluid stress moment induced by the slosh wave 
excitation can be computed from the following relations: 


1 


\FyM,-FM 


_ 1 
^ ^ ^ 






VIII. Methods of Numerical Simulation 

Detailed descriptions of the computational algorithm applicable to 
cryogenic fluid management under microgravity are also given In our earlier 
studies (Hung et al., 1990 a,b,c). In this report, a full-scale GP-B and AXAF 
spacecraft propellant dewar tanks with a radius of 68 cm and a height of 145 cm 
will be used in the numerical simulation. The propellant tank is 80% filled with 
cryogenic liquid helium for the GP-B dewar while the rest of the ullage is filled 
with helium vapor. The temperature of cryogenic helium is 1.8 K. In this study 
the following data were used: liquid helium density - 0,146 g/cm^, helium vapor 
density - 0.00147 g/cm^, fluid pressure - 1.66 x 10^ dyne/cm^, surface tension 
coefficient at the interface between liquid helium and helium vapor - 0.353 
dyne/cm, liquid helium viscosity coefficient - 9.61 x 10*^ cmVs; and contact 
angle - 5°. The initial profiles of the liquid-vapor interface for the rotating 
dewar are determined from computations based on algorithms developed for the 
steady state formulation of microgravity fluid management (Hung et al., 1990 
a,b,c) . 

A staggered grid for the velocity components is used in this computer 
program. The method was developed by Harlow and Welch (1965) for their MAC 
(marker-and-cell) method of studying fluid flows along a free surface. The 
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finite difference method employed in this numerical study was the "Hybrid Scheme" 
developed by Spalding (1972). The formulation for this method is valid for any 
arbitrary interface location between the grid points and is not limited to middle 
point interfaces (Patankar and Spalding, 1972). An algorithm for a semi- implicit 
method (Patankar, 1980) was used as the procedure for modeling the flow field. 
The time step is determined automatically based on the size of the grid points 
and the velocity of flow fields. A detailed description of the computational 
algorithm applicable to microgravity fluid management is illustrated in our 
earlier studies (Hung et al., 1990' a, b,c). Figures 5(A) and 5(B) show the 
distribution of grid points for the dewar tank with probe for GP-B dewar 
container in the radial-axial plane and radial -circumferential plane, 
respectively, in cylindrical coordinates. 

IX. Numerical Simulations of Sloshing Dynamics 
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contract supported by the NASA grant NAGS -129. Results of the research are 
published in the following scientific journals: 

(1) Hung, R. J., Tsao, Y. D., Hong, B. B., and Leslie, F. W. , Axisymmetric 
Bubble Profiles in a Slowly Rotating Helium Dewar Under Low and 
Microgravity Environments, Acta Astronautics . 1£, 411-426, 1989. 

(2) Hung, R. J., Tsao, Y. D. , Hong, B. B., and Leslie, F. W. , Bubble Behaviors 

in a Slowly Rotating Helium Dewar in Gravity Probe-B Spacecraft 

Experiment, J. Spacecraft and Rocke ts. 26, 167-172, 1989. 

(3) Hung, R. J., Tsao, Y. D. , Hong, B. B. , and Leslie, F. W. , Dynamical 
Behavior of Surface Tension on Rotating Fluids in Low and Microgravity 


38 


Environments , International Journal for Mlcrogr avltv Research — an4 
Applications . 11 . 81-95, 1989. 

(4) Hung, R. J., Lee, C. C., and Tsao, Y. D., Fluid Behavior in Microgravity 
Environment, Proceed. National Sci. Council (A) . 14, .21-34, 1990. 

(5) Hung, R. J., Tsao, Y. D., and Leslie, F. W. , Effect of G-Jitters on the 
Stability of Rotating Bubble Under Microgravity Environment, Acta 
Astronautica . 21 . 309-321, 1990. 

(6) Hung, R. J., Tsao, Y. D., Hong, B. B., and Leslie, F. W. , Dynamics of 

Surface Tension in Microgravity Environment, Progress in Aeronautics and 
Astronautics . 127 . 124-150, 1990. 

(7) Hung, R. J. , Superfluid and Normal Fluid Helium II in a Rotating Tank Under 
Low and Microgravity Environments, Proceed. National Sci. Council (A) . 14, 
289-297, 1990. 

(8) Hung, R. J., Lee, C. C. , and Shyu, K. L. , Reorientation of Rotating Fluid 
in Microgravity Environment With and Without Gravity Jitters, J , Spacecraft 
and Rockets . 28 . 71-78, 1991. 

(9) Hung, R. J., Lee, C. C., and Leslie, F. W. , Response of Gravity Level 

Fluctuations on the Gravity Probe-B Spacecraft Propellant System, 2^ 
Propulsion and Power . Z> 556-564, 1991. 

(10) Hung, R. J., Lee, C. C., and Leslie, F. W. , Gravity Jitter Response Slosh 

Wave Excitation on the Fluid in a Rotating Dewar, Advances in Space 

Research . Vol. 11, No. 7, 201-208, 1991. 

(11) Hung, R. J., and Shyu, K. L. , Cryogenic Hydrogen Reorientation and Geyser 
Initiation at Various Liquid- Filled Levels in Microgravity Environment, 
Advances in Space Research . Vol. 11, Mo. 7, 211-226, 1991. 

(12) Hung, R. J., Lee, C. C., and Leslie, F. W. , Gravity- Jitter Effected Slosh 

39 


1 


Waves on the Stability of Rotating Bubble Under Microgravity Environment, 
Advances in Space Research . Vol. 11, No. 7, 209-216, 1991. 

(13) Hung, R. J., and Shyu, K. L. , Cryogenic Liquid Hydrogen Reorientation 
Activated by High Frequency Impulsive Reverse Gravity Acceleration of 
Geyser Initiation, Microeravitv Quarterly . 1, No. 2, 81-92, 1991. 

(14) Hung, R. J., Lee, C. C., and Leslie, F. W. , Slosh Wave Excitation in a 
Partially Filled Rotating Tank Due to Gravity Jitters in a Microgravity 
Environment, Acta Astronautica . 25.. 523-551, 1991. 

(15) Hving, R. J., and Shyu, K. L. , Space-Based Cryogenic Liquid Hydrogen 
Reorientation Activated by Low Frequency Impulsive Thruster of Geyser 
Initiation, Acta Astronautica . 25 . 379-393, 1991. 

(16) Hung, R. J., Shyu, K. L. , and Lee, C. C., Slosh Wave Excitation Associated 
with High Frequency Impulsive Reverse Gravity Acceleration of Geyser 
Initiation, Microeravitv Quarterly , i, No. 3, 125-133, 1991. 

(17) Hung, R. J., Lee, C. C., and Leslie, F. W. , Slosh Wave Excitation of 
Cryogenic Liquid Helium in Gravity Probe-B Rotating Dewar, Advances in 
Crvoeenic Eneineerine . 37, 1281-1290, 1991. 

(18) Hung, R. J., Shyu, K. L. , and Lee, G. C., Cryogenic Liquid Resettlement 
Activated by Impulsive Thrvisc in Microgravity Environment, Advances in 
Cryogenic Engineering . 37. 1313-1326, 1991. 

(19) Hung, R. J., Lee, C. C., and Leslie, F. W. , Gravity Jitters Excited 

Cryogenic Liquid Slosh Waves in Microgravity Environment, Advances is 

Cryogenic Engineering . 37, 1291-1302, 1991. 

(20) Hiing, R. J., Shyu, K. L. , and Lee, C. C. , Slosh Wave Excitation Due to 
Cryogenic Liquid Reorientation in Space-Based Propulsion System, Advances 


40 


In Cryogenic Engineering . 37, 1303-1312. 1991. 


(21) Hung, R. J., and Shyu, K. L. , Constant Reverse Thrust Activated 

Reorientation of Liquid Hydrogen with Geyser Initiation, J . of Spacecraft 
and Rockets . 29 . 279-285, 1992. 

(22) Hung, R. J., and Shyu, K. L. , Fluid Resettlement and Slosh Wave Excitement 

During the Spacecraft Reorientation in Cryogenic Propellant System, 
Proceed. National Scl. Council (4) ■ 131-148, 1992. 

(23) Hung, R. J., Lee, C. C., and Leslie, F. W. , Spacecraft Dynamical 

Distribution of Fluid Stresses Activated by Gravity Jitters Induced Slosh 
Waves, Journal of Guidance. Control and Dynamics , 817-824, 1992. 

(24) Hung, R. J., Shyu, K. L. , Lee, C. C., Liquid Hydrogen Slosh Wave Excited 

by Constant Reverse Gravity Acceleration of Geyser Initiation, jL. 

Spacecraft and Rockets . 29 . 523-528, 1992. 

(25) Hung, R. J. , Shyu, K. L. , and Lee, C. C. , Medium Frequency Impulsive Thrust 
Excited Slosh Waves During Propellant Reorientation with Geyser, J oumal 
of Propulsion and Power . 8, 778-785, 1992. 

(26) Hung, R. J., and Shyu, K. L. , Excitation of Slosh Waves Associated with 
Low Frequency Impulsive Reverse Gravity Acceleration of Geyser Initiation, 
Acta Astronautica . 26 . 425-433, 1992. 

(27) Hung, R. J., and Lee, C. C. , Characteristics and Behaviors of Gravity 
Probe -B Spacecraft Propulsion System, Proceed. Nat ional Science Council 
(A), 16, 239-252, 1992. 

(28) Hung, R. J., and Shyu, K. L. , Suction Dip, Liquid Residual and Liquid 
Propellant Shut-Off in Spacecraft Propulsion System, Proceedings National 
Science Council (4). 1^> 347-364, 1992. 

(29) Hung, R. J., Lee, C. C., and Leslie, F. W. , Similarity Rules in Gravity 


41 


Jitter-Related Spacecraft Liquid Propellant Slosh Waves Excitation, Journal 
of Fluids and Structures . 6, 493-522, 1992, 

(30) Hung, R. J., and Shyu, K. L. , Medium Frequency Impulsive Thrust Activated 
Liquid Hydrogen Reorientation with Geyser, Journal of Propulsion and Power. 
8, 987-994, 1992. 

(31) Htmg, R. J., Lee, C. C., and Leslie, F. W. , Effect of the Baffle on the 
Spacecraft Fluid Propellant Viscous Stress and Moment Fluctuations, 
Transactions of the Japan Society for Aeronautical and Space Sciences . 1^, 
187-207, 1993. 

(32) Hung, R. J., Lee, C. C., and Leslie, F. W. , Dynamics of Gravity Probe-B 
Spacecraft Due to Gravity Jitter Induced Cryogenic Helium Disturbances in 
Rotating Dewar, Advances In Space Research . 13, in press, 1993. 

(33) Hung, R. J., and Shyu, K. L. , Suction Dip, Liquid Residual and Slosh Wave 
Excitation During Liquid Draining in Microgravity, Advances in Space 
Research . 13£Zi, 147-154, 1993. 

(34) Hung, R. J., and Shyu, K. L. , Liquid Settlement, Resettlement, Slosh Wave 
Excitation and Geyser Motion During Reorientation in Microgravity, Advances 
in Space Research . 13 ('7^ . 155-163, 1993. 

(35) Hung, R. J., Lee, C. C., and Leslie, F. W. , Effect of the Baffle on the 
Asymmetric Gravity- Jitter Excited Slosh Waves and Spacecraft Moment and 
Angular Momentum Fluctuations, Journal of Aerospace E ngineering. (United 
Kingdom ) . in press, 1993. 

(36) Hung, R. J., and Pan, H. L. , Asymmetric Slosh Wave Excitation in Liquid- 
Vapor Interface Under Microgravity, Acta Mechanica Slnica . in press, 1993. 

(37) Hung, R. J., and Lee, C. C. , Effect of the Baffle on On-Orbit Spacecraft 


42 


I 

Fluid System Angular Momentum Fluctuations and Action of Viscous Stress 
Moment Exerted on the Container, Proceed. Nation al Science Council (4) . 12. 
in press, 1993. 

(38) Hving, R. J., Pan, H. L. , and Leslie, F. W. , Fluid System Angular Momentum 
and Moment Fluctuations Driven by Gravity Gradient or Gravity Jitter in 
Microgravity, Journal of Flight Sciences a n d Space Research (Zettschrift 
fur Flugwissenschaften und Weltraumf orschung) . in press, 1993. 

(39) Hung, R. J. , and Shyu, K. L. , Liquid Resettlement and Slosh Wave Excitation 
During Fluid Reorientation in Microgravity, Acta Ast ronautica. 22, in 
press, 1993. 

(40) Hung, R. J., Pan, H. L. , and Leslie, F. W. , Gravity Gradient or Gravity 
Jitter Induced Viscous Stress and Moment Fluctuations in Microgravity, 
Fluid Dynamics Research , in press, 1993. 

(41) Hung, R. J., and Shyu, K. L. , Liquid Hydrogen Suction Dip and Slosh Wave 
Excitation Draining Under Normal and Reduced Gravity Environments, 
Transactions of the Japan Society for Aeronautics and Spac e Sciences. 36, 
in press , 1993 . 

(42) Hung. R. J. , and Pan, H. L. , Cryogenic Helium System Angular Momentum and 
Moment Fluctuations Driven by Gravity Gradient Acceleration in 
Microgravity, Advances in Cryogenic Engineering . 38, in press, 1993. 

(43) Hung, R. J., and Lee, C. C., Effect of the Baffle on the Cryogenic Helium 
Container Angular Momentum and Moment Fluctuations Due to Asymmetric 
Gravity Jitter Excited Slosh Waves, Advances n Cryogeni c Engineering. 3£, 
in press, 1993. 

(44) Hung, R. J., and Pan, H. L, , Gravity Gradient Induced Viscous Stress and 


43 


Moment Fluctuations of Cryogenic Helitun System in Microgravity, Advances 
in Cryogenic Engineering . 38 . in press, 1993. 

(45) Hung, R. J., and Shyu, K. L. , Cryogenic Liquid Hydrogen Suction Dip and 
Slosh Wave Excitation During Draining in Microgravity, Advances in 
Cryogenic Engineering. 18. in press, 1993. 

(46) Hung, R. J., and Lee, C. C. , Effect of the Baffle on the Cryogenic Helium 
Container Viscous Stress and Moment Fluctuations In Microgravity, Advances 
in Cryogenic Engineering . 38, in press, 1993, 

(47) Hung, R. J., Lee, C. C., and Leslie, F. W. , Dynamic Characteristics of the 
Partially Filled Rotating Dewar of the Gravity Probe-B Spacecraft, Acta 
Astronautica . 27 . in press, 1993. 

(48) Hung, R. J., and Lee, C. C. , Effect of Baffle on Gravity Gradient 
Acceleration Excited Slosh Waves and Associated Viscous Stress Force 
Activated Spacecraft Dynamic Fluctuations, Aeronautical Journal (Royal 
Aeronautical Society of United Kingdom), in press, 1993. 

(49) Hung, R. J., and Pan H. L. , Differences in Gravity Gradient and Gravity 
Jitter -Excited Slosh Waves in Microgravity, Transactions of the Japan 
Society for Aeronautical and Soace Sciences . 36 . in press, 1993. 

In addition to the publication in the scientific journals, we have also 

published the research results in the national and international scientific 

conferences. The listing of the papers published in the conference proceedings 

are as follows: 

(1) Hung. R. J., Tsao, Y. D. , Hong, B. B. , and Leslie, F. W. , Vibration of 
Microgravity Environment for Rotating Fluids in High and Low Rotating 
Speeds, AIAA Paper 89-0065 . pp. 12, 1989. 

(2) Hung, R. J., Tsao, Y. D., Hong, B. B., and Leslie, F. W. , Dynamics of 


44 


Surface Tension in. Microgravity Environment, Proceedings on Space 
Commercialization , pp. 24, 1989. 

(3) Hung, R. J., Tsao, Y. D. , Hong, B. B., and Leslie, F. W. , Computer Modeling 
of the Dynamics of Surface Tension on Rotating Fluids in Low and 
Microgravity Environments , Finite Element Analysis in Fluids , pp. 211-218, 
1989. 

(4) Hung, R. J., Tsao, Y. D. , Hong, B. B., and Leslie, F. W. , Effect of Gravity 
Level Fluctuations for Rotating Fluids in High and Low Rotating Speeds , 
Finite El ement Analysis in Fluids, pp. 205-210, 1989. 

(5) Hung, R. J., Tsao, Y. D. , and Leslie, F. W. , Effect of G-Jitters on the 
Stability of Rotating Bubble Under Microgravity Environment, Proceedings 
of 40th International Astronautical Congress . lAF Paper 89-4Q9 . pp. 10, 
1989. 

(6) Hung, R. J., Tsao, Y. D. , and Leslie, F. W. , Response of Gravity Level 
Fluctuations on the Gravity Probe-B Spacecraft Propellant System, 
Proceedings of 40th International Astronautical Congress . lAF Paper 89-411 . 
pp. 11, 1989. 

(7) Hung, R. J., Tsao, Y. D. , and Leslie, F. W. , Computer Modeling of the 
Dynamics of Bubble on Rotating Fluids in Low and Microgravity Environments , 
Drops and Bubbles . Ed. by T. G. Wang, pp. 200-207, Published by American 
Institute of Physics, New York, 1989. 

(8) Hung, R. J., Lee, C. C., and Shyu, K. L. , Reorientation of Rotating Fluid 
in Microgravity Environment, AIAA Paper 90-0513 . pp. 12, 1990. 

(9) Hung, R. J., Lee, C. C., Wu, J. L. , and Leslie, F. W. , Gravity Jitters and 
Excitation of Slosh Waves, AIAA Paper 90-0655 . pp. 11, 1990. 

(10) Hung, R. J., Lee, C. C. , and Leslie, F. W. , Slosh Wave Excitation in 


45 


Gravity Probe -B Spacecraft Experiment, Developments In Theoretical and 
Applied Mechanics . H, 197-205, 1990. 

(11) Hung, R. J., Lee, C. C., and Shyu, K. L. , Reorientation of Rotating Fluid 

in Microgravity Environment, Developments in Theoretical and Applied 

Mechanics . 15 . 206-212, 1990. 

(12) Hung, R. J., Lee, C. C. , and Leslie, F. W. , Gravity- Jitters and Excitation 
of Slosh Waves, Developments in Theoretical and Applied Mechanics , 213- 
220, 1990. 

(13) Hung, R. J., Lee, C. C., and Leslie. F. W. , Slosh Wave Excitation and 
Stability of Cryogenic Space Fluid Systems, Worksho p on the 
Commerr-ialization of Space Fluid Management , pp. 194-205, 1990. 

(14) Hung, R. J., Lee, C. C., and Shyu, K. L. , Cryogenic Liquid Propellant 
Reorientations and Fluid Management Used in Space Commercialization, 
Workshop on the Commercialization of Space Fluid Management , pp. 206-215, 
1990. 

(15) Hung, R. J., Lee, C. C., and Shyu, K. L. , Rocket Propellant Reorientation 
and Fluid Management Used in Space Commercialization, Technical and 
Business Exhibition Symposium (TABES) Proceedings , pp. 231-236, 1990. 

(16) Hung, R. J., Lee, C. C., and Leslie, F. W. , Slosh Wave Excitation and 
Stability of Spacecraft Fluid Systems, Technical and Busines s Exhibition 
Symposium (’TABES') Proceedings , pp. 237-244, 1990. 

(17) Hung, R. J., and Shyu, K. L. , Cryogenic Hydrogen Reorientation Activated 
by Constant Reverse Gravity Acceleration of Geyser Initiation, AIAA Paper . 
No. lAF 90-3712 . pp. 10, 1990. 

(18) Hung, R. J., Lee, C. C. , and Leslie, F. W. , Effect of Gravity Jitters on 


46 


the Stability of Rotating Bubble Under Microgravity Environment, tlM 
Paper. No. 91-397 . pp. 11, 1990. 

(19) Hung, R. J. , and Shyu, K. L. , Initiation of Geyser During the Resettlement 
of Cryogenic Liquid in Microgravity Environment, AIAA Pap er. No. 91-0108. 
pp. 11, 1991. 

(20) Hung, R. J., Lee, C. C., and Leslie, F. W. , Gravity- Jitters Excited Slosh 
Waves for Cryogenic Liquid Helium in Rotating Fluid, AIAA Paper. No. 91- 
0488 . pp. 11, 1991. 

(21) Hung, R. J., Lee, C. C. , and Leslie, F. W. , Slosh Wave Excitation in 
Gravity Probe-B Spacecraft Propulsion System, Technical and Business 
Exhibition Symposium (TABES') Proceedings , pp. 54-61, 1991. 

(22) Hung, R. J., and Shyu, K. L. , Cryogenic Liquid Resettlement Activated by 
Impulsive Thrust in Space -Based Propulsion System, Technic al and Business 
Exhibition Symposium (TABES) ProceedinES . pp. 62-70, 1991. 

(23) Hung, R. J., Lee, C. C. , and Leslie, F. W., Gravity Jitters Excited Slosh 
Waves in Rotating Propellant Tank Under Microgravity Environment, Technical 
and Business Exhibition Symposium (TABES) Proceedings , pp. 82-90, 1991. 

(24) Hung, R. J., Shyu, K. L. , and Lee, C. C. , Liquid Propellant Resettlement 
and Slosh Wave Excitation Activated by Impulsive Thrust in Microgravity, 
9th Annual Computational Fluid Dynamics Working Group Meeting, pp. 11, 
1991. 

(25) Hung, R. J., Lee, C. C. , and Leslie, F. W. , Slosh Wave Excitation of 
Cryogenic Liquid Helium in Gravity Probe-B Rotating Dewar, Gryogentd. 
Engineering Conference Proceedinpis . pp. 8, 1991. 

(26) Hung, R, J. , and Shyu, K. L. , Slosh Wave Excitation Due to Cryogenic Liquid 

Reorientation in Space-Based Propulsion System, Cryogenic Engineering 


47 


Conference Proceedings , pp. 9, 1991. 

(27) Hung, R. J., Lee, C. C., and Leslie, F. W. , Gravity Jitters Excited 
Cryogenic Liquid Slosh Waves in Microgravity Environment, Cryogenic 
Engineering Conference Proc eedings, pp. 8, 1991. 

(28) Hung, R. J., and Shyu, K. L. , Cryogenic Liquid Resettlement Activated by 
Impulsive Thrust in Microgravity Environment, Cryogenic Engineering 
Conferenc e Proceedings . pp. 9, 1991. 

(29) Hxing, R. J., Lee, C. C., and Leslie, F. W. , Gravity Probe-B Spacecraft 
Attitude Control Based on the Dynamics of Slosh Wave -Induced Fluid Stress 
Distribution on Rotating Container of Cryogenic Propellant, International 
Astronautical Federation Congress . lAF Paper No. 91-256 . pp. 10, 1991. 

(30) Hung, R. J., Shyu, K. L. , and Lee, C. C., Liquid Resettlement and Slosh 
Wave Excitation During Fluid Reorientation in Microgravity, Interna tional 
Astronautical Federation Congress . lAF Paper No. 91-392 . pp. 11, 1991. 

(31) Hung, R. J., Lee, C. C., and Leslie, F. W. , Longitudinal and Transverse 
Modes of Slosh Wave Excitation in Rotating Dewar Associated with Gravity 
Jitters, International Astronautical Federation Congress . lAF Paper No. 91- 

f 

404 . pp. 11, 1991. 

(32) Hung, R. J., Pan, H. L. , Lee, C. C., and Leslie, F. W. , Effect of 
Asymmetric Gravity Jitter Excited Slosh Waves at Liquid-Vapor Interface 
Under Microgravity, AIAA Paper . No. 92-0247 . pp. 11, 1992. 

(33) Hung, R. J., and Shyu, K. L. , Draining Characteristics of Cryogenic Liquid 
in a Low Gravity Environment, AIAA Paper . No. 92-0819 . pp. 10, 1992. 

(34) Hung, R. J., Lee, C. ~C., and Leslie, F. W. , Fluid Stress Disturbance 
Activated by Gravity Jitter Induced Slosh Waves in Microgravity, AIAA 
Paper . No. 92-0880 . pp. 10, 1992. 


48 


(35) Hung, R. J. , and Shyu, K. L. , Numerical Slmulaclon of Spacecraft Propellant 
Tank Draining Characteristics, Developments in Theoretical and Applied 
Mechanics . 12, II . 9 . ll-II . 9 . 23. Nashville, Tennessee, 1992. 

(36) Hung, R. J., Lee, C. C., and Leslie, F. W. , Similarity Rules in Testing of 
Microgravity Fluid Behaviors, Developments in Theoretical and Applied 
Mechanics . 17 . II . 9 . l-II . 9 . 10, Nashville, Tennessee, 1992. 

(37) Hung, R. J., Pan, H. L. , and Leslie, F. W. , Three-Dimensional Spacecraft 
Angular Momentum Fluctuations Due to Gravity Jitter Excited Slosh Waves in 
Liquid Propellant System, International Astronautical Federation Congress . 
lAF Paper No. 92-0025 . pp. 11. 1992. 

(38) Hung, R. J., Lee, C. C. , Leslie, F. W. , Gravity Jitter Induced Slosh Waves 
in Liquid Propellant and Disturbances on Spacecraft Angular Momentum, 
International Astronautical Federation Congress , lAF Paper No. 92-0667 , pp . 

10, 1992. 

(39) Hung, R. J. , Shyu, K. L. , Liquid Draining Shut-Off Induced Geyser and Slosh 
Wave Excitation on Suction Dip During Draining in Microgravity, 
International Astronautical Federation Congress . lAF Paper No. 92-0921 . pp. 

11, 1992. 

(40) Hung, R. J., Pan, H. L. , and Leslie, F. W. , Asymmetric Gravity Jitter 
Excited Slosh Waves at Liquid-Vapor-Solid Interface Under Microgravity, 
International Astronautical Federation Congress . lAF Paper No. 92-0920 . 
pp. 10, 1992. 

(41) Hung, R. J., Lee, C. C., and Leslie, F. W. , Dynamics of Gravity Probe-B 
Spacecraft Due to Gravity Jitter Induced Cryogenic Helium Disturbances in 
Rotating Dewar, COSPAR . 29, pp. 14, 1992. 

(42) Hung, R. J., and Shyu, K. L. , Suction Dip, Liquid Residual and Slosh Wave 


49 



Excitation During Liquid Draining in Microgravity, COSPAR . 2£, pp. 15, 
1992. 

(43) Hung, R. J., and Shyu, K. L. , Liquid Settlement, Resettlement, Slosh Wave 
Excitation and Geyser Motion During Reorientation in Microgravity, COSPAR. 
21 . pp. 1^. 1992. 

(44) Hung, R. J., Lee, C. C. , and Leslie, F. W. , Effect of the Baffle on 
Asymmetric Gravity Jitter Excited Slosh Waves in Microgravity, AIAA Paper . 
No. 93-0913 . pp. 11, 1993. 

(45) Hung, R. J., and Pan, H. L. , Cryogenic Helitun System Angular Momentum and 
Moment Fluctuations Driven by • Gravity Gradient Acceleration in 
Microgravity, Cryogenic Engineering Conference Proceedings , in press, 1993. 

(46) Hung, R. J., and Lee, C. C. , Effect of the Baffle on the Cryogenic Helium 
Container Viscous Stress and Moment Fluctuations in Microgravity , Cryogenic 
Engineering Conference Proceedings , in press, 1993. 

(47) Hung, R. J., and Pan, H. L. , Gravity Gradient Induced Viscous Stress and 
Moment Fluctuations of Cryogenic Helitun System in Microgravity, Cryogenic 
Engineering Conference Proceedings . in press, 1993. 

(48) Hung, R. J., and Shyu, K. L. , Cryogenic Liquid Hydrogen Suction Dip and 
Slosh Wave Excitation During Draining in Microgravity, Cryogenic 
Engineering Conference Proceedings , in press, 1993. 

(49) Hung, R. J., and Lee, C. C. , Effect of the Baffle on the Cryogenic Helium 
Container Angular Momentum and Moment Fluctuations Due to Asymmetric 
Gravity Jitter Excited Slosh Waves, Cryogenic Engineering Conference 
Proceedings . in press, 1993. 

(50) Hung, R. J., and Leslie, F. W. , Characteristics of Gravity Probe-B 
Spacecraft Cryogenic Propellant Management, Proceedings of First Pacific 


50 


International Conference on Aerospace Sc ience and Technology. In press, 
1993. 

(51) Hung, R. J. , and Pan, H. L. , Acting of Gravity Gradient and Gravity Jitter 
Accelerations on the On-Orbit Orbital Spacecraft, Proceed ings of First 
Pacific International Conference on Aerospace Science and Technology, in 
press, 1993. 

(52) Hung, R. J., and Shyu, K. L. , Liquid Propellant Residual and Shut-Off in 
Spacecraft Propulsion System, Proceedings of First Pacific I nternational 
Conference on Aerospace Science and Technology , in press, 1993. 

(53) Hung, R. J., and Lee, C. C., Effect of the Baffle on On-Orbit Spacecraft 
Angular Momentum and Viscous Moment Exerted on Propellant System, 
Proceedings of First Pacific International Conference on Aerospace Science 
and Technology , in press, 1993. 

(54) Hung, R. J., and Shyu, K. L. , Vapor ingestion. Liquid Residual and Slosh 
Wave Excitation During Two-Phase Flow Drainage and Shut-Off in 
Microgravity, AIAA Paper . No. 94-0703 . pp. 11, 1994. 

(55) Hung, R. J., and Pan, H. L. , Liquid-Vapor Interface Fluctuations Induced 
by Gravity Gradient and Jitter Accelerations in Microgravity, AIAA Paper . 
No- 94-0439 . pp. 12, 1994. 

Acknowledgement 

The authors appreciate the support received from the National Aeronautics 
and Space Administration through the NASA Grant NAG8-129. They would like to 
express their gratitude to Richard A. Potter and Fred W. Leslie of NASA/Marshall 
Space Flight Center for the stimulating discussions during the course of the 
present study. 


51 


References 


Avduyevsky V. S. (editor), Scientific Foundations of Space Manufacturing, 

MIR, Moscow, USSR 1984. 

Everitt, F. Ed., Testing Einstein With Orbiting Gyroscopes - Gravity P robe-B. pp. 
28, Stanford University, Palo Alto, CA, 1990. 

Forward, R. L. , Flattening Space-Time Near the Earth, Physical Review . Series 2, 
Vol. 26, pp. 735-744, 1982. 

Hoare, F. E., Jackson, L. C., and Kurti, N. , Experimental Cryogenics : Liquid 
HeTfum TT . Butterworths , London, U.K. (1961). 

Hung, R. J., Superfluid and Normal Fluid Helium II in a Rotating Tank Under Low 
and Microgravity Environments, Proceedings of National Science Council . 
Series (A) . Vol. 14, pp. 289-29, 1990. 

Hving, R. J., and Lee, C. C., Characteristics and Behaviors of Gravity Probe-B 
Spacecraft Propulsion System, Proceed. National Science Council (A) . 16, 
pp. 339-352, 1992. 

Hung, R. J. , Lee, C. C. , and Leslie, F. W. , Dynamics of Gravity Probe-B 
Spacecraft Due to Gravity Jitter Induced Cryogenic Helium Disturbances in 
Rotating Dewar, Advances in Space Research . Vol. 13, in press, 1993c. 

Hung, R. J., Lee, C. C. , and Leslie, F. W. , Effect of G-Jitters on the Stability 
of Rotating Bubble Under Microgravity Environment, Acta Astronautica . Vol. 
21, pp. 309-321, 1990a. 

Hung, R. J., Lee, C. C. , and Leslie, F. W. , Gravity- Jitter Response Slosh Waves 
Excitation on the Fluid Rotating Dewar, Advances in Space Research . Vol. 

- 11(7), pp. 201-208, 1991d. 

Hung, R. J., Lee, C. C. , and Leslie, F. W. , Gravity- Jitter Effected Slosh Waves 
on the Stability of Rotating Bubble Under Microgravity Environment, 


52 


Advances tn Space Research . Vol. 11(7), pp. 209-216, 1991e. 

Hung, R. J., Lee, C. C. , and Leslie, F. W. , Slosh Wave Excitation in a Partially 
Filled Rotating Tank Due to Gravity Jitters in a Microgravity Environment, 
Acta Astronautica . Vol. 25, pp. 523-551, 1991f. 

Hung, R. J., Lee, C. C. , and Leslie, F. W. , Response of Gravity Level 
Fluctuations on the Gravity Probe-B Spacecraft Propellant System, Journal 
of Propulsion and Power . Vol. 7, pp. 556-564, 1991g. 

Hung, R. J., Lee, C. C. , and Leslie, Slosh Wave Excitation of Cryogenic Liquid 
Helium in Gravity Probe-B Rotating Dewar, Advances in Cryogenic 
Engineering . 37 . 1281-1290, 1991J . 

Hung, R. J., Lee, C. C., and Leslie, F. W. , Gravity Jitter Excited Cryogenic 
Liquid Slosh Waves in Microgravity Environment, Advances in Cryogenic 
Engineering . 37, 1291-1302, 1991c. 

Hung, R. J. , Lee, C. C. , and Leslie, F. W. , Spacecraft Dynamical Distribution 
of Fluid Stresses Activated by Gravity Jitter Induced Slosh Waves, Journal 
of Guidance . Control and Dynamics . Vol. 15, pp. 817-824, 1992c. 

Hung, R. J., Lee, C. C., and Leslie, F. W. , Similarity Rules in Gravity Jitter - 
Related Spacecraft Liquid Propellant Slosh Waves Excitation, Journal of 
Fluids and Structures . Vol. 6, pp. 493-522, 1992d. 

Hung, R. J., Lee, C. C. , and Leslie, F. W. , Effect of the Baffle on the 
Spacecraft Fluid Propellant Viscous Stress and Moment Fluctxiations , 
Transaction of the Japan Society for Aeronautical and Space Sciences . 
187-207, 1993a. 

Hung, R. J. , Lee, C.C. , and Leslie, F. W. , Effect of the Baffle on the Asymmetric 
Gravity- Jitter Excited Slosh Waves and Spacecraft Moment and Angular 
Momentum Fluctuations , Journal of Aerospace Engineering (United Kingdom) , 


53 


In press, 1993b. 

Hung, R. J., Lee, C. C. , and Shyu, K. L. , Reorientation of Rotating Fluid in 
Microgravity Environment With and Without Gravity Jitters, Journal of 
Spacecraft and Rockets . Vol. 28, pp. 71-79, 1991h. 

Hung, R. J., and Leslie, F. W. , Bubble Shapes in a Liquid-Filled Rotating 
Container Under Low Gravity, Journal of Spacecraft and Rockets. Vol. 25, 
pp. 70-74, 1988. 

Hung, R. J., Pan, H. L. , and Leslie, F. W. Gravity Gradient or Gravity Jitter 
Induced Viscous Stress and Moment Fluctuations in Microgravity, Fluid 
hvTianilcs Research , in press 1993d. 

Hung, R. J., and Pan, H. L. , Asymmetric Slosh Wave Excitation in Liquid-Vapor 
Interface Under Microgravity, Acta Mechanica Sinica . in press, 1993. 

Hung, R. J., and Shyu, K. L. , Cryogenic Hydrogen Reorientation and Geyser 
Initiation at Various Liquid- Filled Levels in Microgravity, Advances in 
Space Research . Vol. 11(7), pp. 217-226, 1991a. 

Hung, R. J., and Shyu, K. L. , Cryogenic Liquid Hydrogen Reorientation Activated 
by High Frequency Impulsive Reverse Gravity Acceleration of Geyser 
Initiation, Microgravitv Quarterly . Vol. 1(2), pp. 81-92, 1991b. 

Hung, R. J. , and Shyu, K. L. , Space-Based Cryogenic Liquid Hydrogen Reorientation 
Activated by Low Frequency Impulsive Reverse Gravity Thruster of Geyser 
Initiation, Acta Astronautics . Vol. 25, pp. 709-719, 1991c. 

Htmg, R. J., and Shyu, K. L. , Constant Reverse Thrust Activated Reorientation 
of Liquid Hydrogen with Geyser Initiation, Journal of S pacecraft and 
Rockets ■ Vol. 29, pp. 279-285, 1992a. 

Hung, R. J., and Shyu, K. L. , Excitation of Slosh Waves Associated with Low 
Frequency Impulsive Reverse Gravity Acceleration of Geyser Initiation, 


54 


Astronautica . Vol. 26, pp. 425-433, 1992b. 

Hung, R. J., and Shyu, K. L. , Medium Frequency Impulsive Thrust Activated 

Liquid Hydrogen Reorientation with Geyser, Journal o f Propulsion and Power. 
Vol. 8, pp. 987-994, 1992c. 

Hung, R. J., and Shyu, K. ., Suction Dip, Liquid Residual and slosh Wave 

Excitation during Liquid Draining in Microgravity, Advances Ip Spage 

Research . 13(7) . 147-154, 1993a. 

Hung, R. J., and Shyu, K. L. , Liquid Settlement, Resettlement, Slosh Wave 
Excitation and Geyser Motion During Reorientation in Microgravity , Advances 
in Space Research . 13(7) . 155-163, 1993b. 

Hung, R. J., and Shyu, K. L. , Liquid Resettlement and Slosh Wave Excitation 
During Fluid Reorientation in Microgravity, Acta Ast ronautica. 27, in 
press, 1993c. 

Hung, R. J., Shyu, K. L. , and Lee, C. C., Slosh Wave Excitation Associated with 
High Frequency Impulsive Reverse Gravity Acceleration of Geyser Initiation, 
Microeravitv Quarterly . Vol. 1(3), pp. 125-133, 19911. 

Hung, R. J., Shyu, K. L. , and Lee, C. C. , Cryogenic Liquid Resettlement Activated 
by Impulsive Thrust in Microgravity Environment, Advances in Cryogenic 
Engineering . 37, 1313-1326, 1991a. 

Hung, R. J., Shyu, K. L. , and Lee, C. C., Slosh Wave Excitation Due to Cryogenic 

Liquid Reorientation in Space-Based Propulsion System, Advances in 

Cryogenic Engineering . 37, 1303-1312, 1991b. 

Hung, R. J., Shyu, K. L. , and Lee, C. C., Medium Frequency Impulsive Thrust 

Excited Slosh Waves During Propellant Reorientation with Geyser , J ournal 
of Propulsion and Power . Vol. 8, pp. 778-785, 1992a. 

Hung, R. J., Shyu, K. L. , and Lee, C. C. , Liquid Hydrogen Slosh Wave Excited by 


55 


Constanc Reverse Gravity Acceleration of Geyser Initiation, Journal of 
Spacecraft and Rockets . Vol. 29, pp. 523-528, 1992b. 

Hung, R. J., Tsao, Y. D., Hong, B. B. , and Leslie, F. W. , Time Dependent 

Dynamical Behavior of Surface Tension on Rotating Fluids Under Microgravity 
Environment, Advances in Space Research . Vol. 8(12), pp. 205-213, 1989d. 

Hung, R. J., Tsao, Y. D., Hong, B. B., and Leslie, F. W. , Dynamics of Surface 

Tension in Microgravity Environment, Progress in Aeronautics — and 
Astronautics . Vol. 127, pp. 124-150, 1990b. 

Hung, R. J., Tsao, Y. D., Hong, B. B., and Leslie, F. W. , Dynamical Behavior of 
Surface Tension on Rotating Fluids in Low and Microgravity Environments, 
International Journal for Microeravitv Research and Applications . Vol. 11, 
pp. 81-95, 1989a. 

Hung, R. J., Tsao, Y. D. , Hong, B. B. , and Leslie, F. W. , Axisymmetric Bubble 
Profiles in a Slowly Rotating Helium Dewar Under Low and Microgravity 
Environments, Acta Astronautica . Vol. 19, pp. 411-426, 1989b. 

Hung, R. J., Tsao, Y. D. , Hong, B. B., and Leslie, F. W. , Bubble Behaviors 

in a Slowly Rotating Helium Dewar in Gravity Probe-B Spacecraft Experiment, 
Journal of Spacecraft and Rockets . Vol. 26, pp. 167-172, 1989c. 

Kamotani, Y. , Prasad, A., and Oastrach, S., Thermal Convections in an Enclosure 
Due to Vibrations Aboard a Spacecraft, AIAA Journal . Vol. 19, pp. 511-516, 
1981. 

Landau, L. D. , and Lifshitz, Fluid Mechanics . Pergamon Press, London, pp. 1-656, 
1959. 

Leslie, F. W. , Measurements of Rotating Bubble Shapes in a Low Gravity 
Environment, Journal of Fluid Mechanics . Vol, 161, pp. 269-275, 1985. 

Mason, P. , Collins, D. , Petrac, D., Yang, L. , Edeskuty, F,, Schuch, A., and 


56 


Williamson, K. , The Behavior of Superfluid Helium in Zero Gravity, 
Proceedings 7th International Cryogenic Eneln eerinE Conferenceg. Surrey, 
England, Science and Technology Press, 1978, 

Misner, C. W. , Thome, K. S., and Wheeler, J. A., Gravitation. W. H. Freeman 
Co., San Francisco, CA, pp. 1-1279, 1973. 

Stanford Relativity Gyroscope Experiment (NASA Gravity Probe-B), Proceedings of 
Society of Photo-Ootical Instrumentation Engineers . £19, pp. 1-165, Society 
of Photo-Optical Instrumentation Engineers, Bellingham, WA, 1986. 

Weinberg, S., Gravitation and Cosmology - Principles and Appli cations of General 
Relativity , pp. 657, John Wiley and Sons, New York, 1972. 

Wilks, J., The Properties of Liquid and Solid Helitim . Clarendon Press, Oxford, 
U.K., 1967. 

Wilkinson, D. T. , Bender, P. I., Eardley, D. M. , Gaisser, T. K. , Hartle, J, B., 
Israel, M. H. , Jones, L. W. , Partridge, R. B. , Schramm, D. N., Shapiro, I. 
I., Vessort, R. F. C. , and Wagoner, R. V., Gravitation, Cosmology and 
Cosmic Ray Physics, Physics Today . 19, 43-46, 1986. 


57 



Figure Captions 


Figure 1 


Figure 2 


Figure 3 


Figure 4 


Figure 5 


GP-B spacecraft coordinate systems with azimuth angle ^ from 
spacecraft mass center to the center of the Earth. Coordinate (x" , 
y" , z") for slew motion and coordinate(x, y, z) for fluid mechanics 
computation. 

Geometry of the GP-B dewar container with the coordinate system 
perpendicular and tangential to the container wall. (A) Geometry in 
r-z plane and (B) Geometry in r-tf plane. 

Geometry of the GP-B dewar container with baffle-boards and their 
locations. (A) Geometry in r-z plane and (B) Geometry in r-^ plane. 
The GP-B module showing main elements of liquid helium dewar, probe 
and baffle-boards. 

Distribution of grid points in the (A) Radial-axial plane, and (B) 
the radial -circtimferential plane of the cylindrical coordinates for 
the GP-B dewar tank. 
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